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C. D. U. 624.313 


LES LIAISONS INTERNES DES MACHINES 
ELECTRIQUES A COMMUTATION 


par MANUEL S. GARRIDO 


Boursier de II, A. €,, Lisbonne 
Lab. Electrotechnique, 
E.N,S. E.M,, Nancy 


RESUMO SYNOPSIS 

Mostra-se como estabelecer as equações gerais das má- It is shown how to establish the general equations of 
quinas eléctricas de comutação, tendo em consideração as commutator machines, including in the analysis the commu- 
secções em comutação e como estabelecer um modelo equiva- tating coils and how to establish an equivalent model where 
lente à máquina real no qual os efeitos das correntes de the effects of commutating currents are negligeable, 


comutação sejam desprezáveis. 


|. INTRODUCTION 


On a montré dans des publications précédentes que la loi d'induction sous la forme e = — dy/dt 
n'est pas directement applicable aux machines électriques à commutation en fonction uniquement 

des variables et des parametres externes de la machine [1] [2] |3]. 

En fait, si Ion veut appliquer la loi d'induction sous forme e = — d/dt il faut introduire dans 
[analyse du systême les sections en commutation. On a déjá montré [3] comment cette méthode 
permet d'obtenir les équations d'une machine à courant continu simple, non saturée. 

On a alors vérifié que les équations de fonctionnement des machines électriques à commu- 
tation sont des équations de valeurs moyennes prises sur l'intervalle de commutation. 

Ceci nous a permis d'établir un modêle équivalent à la machine à commutation, capable de 
représenter ce qui se passe en valeur moyenne sur lintervalle de commutation dans la machine 
réelle. Les équations de fonctionnement de ce modêle équivalent seront, alors, les équations cherchées 
pour la machine réelle. On a ainsi pu établir les équations générales des machines à commutation, 
saturées ou non, avec un nombre quélconque de circuits électriques [3]. 

Le modéle équivalent utilisé est de même nature que la roue de Barlow et pour I'établisse- 
ment de ses équations de fonctionnement, on ne peut plus utiliser la loi d'induction sous la forme 
e=— dy'dt. D'un point de vue dynamique, cela signifie que les équations de Lagrange ne sont pas 
applicables. En réalité, on a vérifié qu'il faut tenir compte de certaines relations entre les variables du 
modéle équivalent, relations qui nous permettent d'obtenir équations dynamiques du systéme. 

Dans le présent article, nous nous proposons d'analyser de plus prês le passage de la machine 
réelle à son modele équivalent. 

Ainsi, nous commencerons par montrer comment établir les équations générales des machines 
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électriques à commutation en introduisant les sections en commutation, ce qui revient à travailler 
dans le domaine oú les équations de Lagrange sont applicables. 

Aprês, nous allons vérifier comment le passage à des valeurs moyennes conduit directement 
au modéle éguivalent et aux relations entre variables dont on a besoin pour établir ses équations 
dynamiques 

On supposera les machines dépourvues d'hysteresis et on négligera les courants de Foucault. 
On négligera, également, I'énergie électrostatique. 

Dans lanalyse qui va suivre on admettra I'hypothese d'avoir plusieures machines en série 
ce qui ne complique pas I'exposé. 


IH. ETUDE DE LA COMMUTATION POUR LA THEORIE DYNAMIQUE DES MACHINES 
ELECTRIQUES 


Désignons par a, b... les circuits électriques independants, d'un systêéme électromécanique. 
Ces circuits pourront être à commutation ou non et un même circuit pourra être distribué sur plu- 
sieurs armatures. 

Désignons par O, les coordonnées généralisées des armatures et par pau les coordonnées de 
position des «balais» des circuits a mesurées selon les coordonnées (O, [3]. Une coordonnée au est 
définie seulement quand le circuit a sétablit dans l'armature concernant la coordonnée O. 

On supposera partout que les machines à I'étude sont bipolaires, symétriques et leurs enrou- 
lements diamétraux, à commutation simple. 

Désignons par iau le courant de commutation d'un circuit a associé à la coordonnée (,, les 
deux sections en commutation étant supposées parfaitement couplées entre elles. 

L'équation électrique instantanée d'un circuit a sera: 


d' ' d dT'E 
(1) va = Ra ia + Ya = Raia + E 
dt dt dia 
ou T'g désigne la coénergie magnétique du systême. 
Cette équation peut être développée comme suit : 
y )2 “E a , )2 é j , 4 Jd T'p 1 
(2) n=Rin+AS + õ. E 0 Te Mo A E Te di 


u 
u diadOy D diadip dt bb dia Odipu dt 


Une machine à commutation est caractérsée par le fait que les phénomeênes se répetent péri- 
odiquement, avec une période qu'on désignera par Atc. Les caractéristiques de la machine corres- 
pondent à son comportement movyen sur lintervalle Atc. 

Quand il y a un seul circuit à commutation, la période Atc correspond à un seul intervalle de 
commutation. Lorsqu'il y a plusieurs circuits à commutation, la période At. sera le plus petit com- 
mun multiple des différents intervalles de commutation Átau. 

La valeur moyenne de Iéquation (2) sur Vintervalle Átç s'écrit 


to + atc 


t te am é 
(3) + | vadtm o | (Raia +3 oTe. Oy + 


Ate to Àte to u dia dO 
sz E 9º TE dip lh Y o? PÉ 2) 


— 


b dia dip dt b,u dia dibu dt 


Supposons que nous ayons des machines avec un nombre três élevé de sections dans les en- 
roulements à commutation. Ce fait nous permet plusieurs simplifications. 

On peut alors admettre que les valeurs moyennes, sur I'intervalle tc, des courants électriques 
ia et des vitesses mécaniques O, ne diffêrent pratiquemeut pas de leurs valeurs instantanées. 
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En plus, la contribution des courants de commutation iau pour la coénergie magnétique totale 
sera négligeable. 
Alors, I'équation (3) pourra s'écrire 


(4) va=Raia + 2- é l 
o O 4 dia 00, b diadip dt bu diadipu Bte u 


J? Te. o, n d* Te  dip +s o Pe 1 [o Fate dibu 
to dt 


Puique lintervalle Atc est un multiple des intervalles Átau, on peut affirmer que: 


(5) 


1 to + ate di i 

[ bu Ra 598 b 

Ate “ to dt Atpe 

lorsqu'on considere des enroulements avec deux voies en paralléle. 
L'équation (4) peut alors s'écrire:. 


0º T' : ge Te di TR 
(6) Ema A, As E sã EA + 
u dia dO b dia dipb dt b,u d ia dibu A tbu 


Une transformation immédiate nous permet maintenant d'écrire I'équation (6) sous la forme: 


7 ( dT'E 


: 9º Te dip dT'E ib 1 dT'E 
7 va=Raia + 2 >—— d — e ) s 
7) i capas b diadip dt bu dia N dO, : dibu Atbu u ÀÁtau diau 


Pour analyser cette éguation, fixons notre attention sur un circuit à commutation b. 

Les conducteurs de ce circuit sont entrainés par les armatures avec des vitesses O, et les temps 
de commutation du circuit b sont Átpu. 

Supposons maintenant que pendant un certain intervalle de commutation Atpy tout le systeme 


électromécanique soit immobilisé à exception du circuit b qui conservera sa vitess Opy = O, selon 
la coordonnée O, . Pour cela, ses conducteurs seront dégagés des liaisons mécaniques avec l'armature 
concernant la coordonnée Oy. 

En plus, on admettra que touts les courants sont maintenus constants, à l'exception du courant 
ibv lequel, par suite des conditions de liaison du systême sera obligé de s'inverser pendant Vinter- 
valle Atpy. Cette inversion se réalise, naturellement, pendant lintervalle Atpy, ou alors forcément 
à la fin de cet intervalle. 

Le résultat de linversion du courant iby est que la configuration du systême des courants 
ip et ipy rest la même aux instants to et to + 4tby. Alors, aux instants to et to + Átpy la distribution 
de champ magnétique dans le systême et sa coénergie magnétique seront les mêmes. 

Cette conclusion nous permet d'affirmer que: 


* to + Atby O TE 4 od T'e; dipy 
a ( ELE Qd =) dt=o0 
) ) to o) Opy és 7) iby dt 


D'aprês les remarques déja faites, (8) peut s'écrire sous la forme: 


"4 DR. TI 


d py diby Atpy 


(9) 


Supposons maintenant que les balais du circuit b sont mobiles et que pendant lintervalle de 
temps À tpy, on a conservé leur vitésse opy. 
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Dans ce cas, la coénergie magnétique du systême aura changé entre les instants to et to + À thy 
d'une quantité 


(10) ATE = — 


Alors, dans le cas général des balais tournants, I'équation (9) prendra la forme: 


Te. dTE à o dTE: 


(11) Om — E = E, 
dO, ni dipy A tpy dOpy cá 


Les résultats obtenus pour le circuit b seront valables pour tous les circuits dans le mouvement 
réel du systême, ce qui nous permet d'écrire: 


(12) 3 


ss E d Te ib )= o Te . “TE: 
b,u 


d Oy li d ibu Â Ebu 


ayant introduit une variable O,, qui tient compte des variations de reluctance magnétique associées 
au mouvemente de l'armature selon la coordonée Ou. 

L'équation (12) nous permet de modifier l'avant-dernier terme d'équation (7). Le dernier terme 
de cette équation peut également être écrit sous une forme plus commode, obtenue à partir de (11): 


(13) 4 Te: Pau 2 Te 
u Átau Odiau u da dOau 


ayant introduit la désignation 1, pour représenter les vitesses (0, = cau)- 
Ainsi, I'équation (7) prendra la forme finale: 
9º T'e 


"Te dip Ta . 1 
(14 v=Rits>—E4+s—o +tI— + q 
ú pi b dia dip dt u diad Ouy » bu Oia dOpu di u ia dO 


Les équations (14) posseêdent, par rapport aux équations (7) Iavantage de ne contenir que des 
variables externes du systême. Les courants de commutation ont été éliminés et ils n'interviendront 
plus dans V'étude du mouvement de la machine, puisque leur contribution pour la coénergie magné- 
tique T'p est négligeable. 

Les équations ainsi obtenues pour le cas général des machines à commutation à balais tour- 
nants contiennent comme cas particulier celui des machines sans commutation. Celles-ci correspon- 


dent au cas oú toutes les vitesses nau sont nulles et les vitesses O et Pau égales. 
Ainsi, l'équation électrique d'un circuit sans commutation sera: 
9º T'e dib Te . 9º T'g 


(15) Va== Ra ia + 3 —— à A qo 
' = a dia dip dt + u dia dOuu ud EA dia d Opy e 


qu'on pourra écrire plus simplement : 


2Te dir, s PTe . 
(16) n=Riti-—— 2+2-— 4=Rito— E 
“E O digdis dt * waão, a dt dia 


Les équations électriques des machines à balais fixes sont obtenues par contre, quand toutes 
les vitesses cau sont nulles, les Oy et nau étant alors égaux. 
Les équations mécaniques des machines à commutation sont établies sans difficulté. En fait, 
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les forces électromagnétiques ne dépendent que de la configuration du systême et des valeurs des 
courants. Ainsi, les forces électromagnétiques seront toujours données par : 


dT' 
17 r si 
Ara Da dO 


que les circuits électriques soient à commutation ou non. 
Si on veut séparer les forces électrodynamiques exercées sur les circuits électriques, des forces 
dues aux variations de reluctance, on peut écrire: 


oT'e d TE 
(18) ( o) à Ouu cá à Ou 


/ em 
On remarquera que la contribution des sections en commutation à la force (['u)em est négligeable. 
Exemple 


Considérons le cas d'une machine à courant continu classique, avec un seul circuit inducteur, 
désigné par 1 et seul circuit induit désigné par 2. Soit O Iangle de position du rotor. 
L'équation de linducteur sera d'aprês I'équation (14): 


du du | du dr 


du dt dis dt 


vu=Rysiu+ 


L'équation d'induit prend la forme: 


Wa j ) ya di 6) 9d TE 
va== Raiz Bo du. + a dp é 
du dt dia dt ii dO 
Enfin, le couple électromagnétique est: 
r Te 
a" d 6) 


Si on considêre le régime permanent, on verifie tout de suite que, une fois déduites les pertes 
Joule, toute I'énergie électrique reçue par Winduit est convertie en énergie mécanique. 
Considérons l'expression de la coénergie magnétique : 


[Eta (iu, io, O) dio 


i 
“1,0 


Te= [0º bd (is, o, O) di + 
o 


o, 


ow (i,0,0) 


Supposons que dV/d O est indépendant de iz. Alors, puisque 0 


==0 avec un rotor lisse, 


odT'e + ka o ya 


dO 9 O 


Ce résultat nous permet d'écrire l'éguation d'induit, en régime permanent 


J Va Ê 
vo == Ra is (=) 
i + d €) 


On arrive ainsi à une expression três répandue de la f. e. m. d'induit d'une machine à courant continu. 

Pourtant cette équation n'est généralement correcte que pour la machine à vide (iz==0), ou 
pour une machine non saturée avec les balais calés sur un axe de symétrie, car alors d bs/dO est 
indépendant de iz. 
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HI. LIAISONS INTERNES DANS LES MACHINES A COMMUTATION 


Quand on considêre le comportement instantané d'une machine à commutation, à l'intérieur 
de Vintervalle de commutation, les équations de Lagrange sont applicables et conduisent aux égua- 
tions électriques (2) et aux équations mécaniques (17). Les f.e.ms. et les forces électrodynamiques 
sont alors calculées comme des forces d'inertie généralisées. 

Pourtant, les équations qui nous intéressent dans les machines à commutation sont, en général, 
les équations des valeurs moyennes prises sur l'intervalle de commutation. Ainsi, le point de vue 
qui nous intéresse pour analyser le comportement des machines à commutation est celui d'un obser- 
vateur insensible à ce qui se passe à lintérieur de lintervalle de commutation. Cet observateur 
tiendra compte seulement du comportement moyen de la machine et ainsi, il doit être capable 
d'établir directement les équations de fonctionnement (14), (17). 

Il nous faut alors savoir comment cet observateur décrit les phénomeênes et comment il établit 
les équations de fonctionnemet. 

La premiére conclusion importante concerne les dérivés des fonctions d'énergie par rapport 
aux coordonées mécaniques. 

Ainsi, considérons un seul enroulement à commutation porté par une armature lisse: p.e., 
linduit d'une machine à courant continu classique. 

Supposons qu'on fasse tourner larmature d'une lame de collecteur ou d'un nombre entier de 
lames de collecteur, les courants extérieurs de la machine étant maintenus constants. 

L'énergie magnétique de la machine a souffert des petites variations instantanées, auxquelles 
observateur des valeurs moyennes n'est pas sensible, et reprend à la fin sa valeur initiale. 

Alors, I'observateur sera conduit a affirmer que: 

dTe O Te 


(19) q O o Ê d E 


+ . 4 Fed 4 ) * . . . 4 + 
[l est évident que les dérivés e ont des significations différentes pour I'observateur des 
dt 


valeurs instantanées et pour l'observateur des valeurs moyennes. Pour le premier, il s'agit de dérivés 
prises avec les courants iau constants, alors que pour le deuxiême les dérivées incluent I'effet du 
changement des courants iau qui est imposé par les liaisons de la machine. 

Si la machine est linéaire et qu'on puisse définir des paramétres, les conditions (19) signifient 
que pour lobservateur des valeurs moyennes, ces parameétres sont indépendants des coordonées mé- 
caniques €), C'est-a-dire, il prendra comme paramêtre de la machine les valeurs moyennes sur Vin- 
tervalle de commutation des parametres mesurés par l'observateur des valeurs instantanées. 

Supposons maintenant que,en plus d'un déplacement de l'armature, I'observateur des valeurs 
moyennes donne aussi un déplacement aux balais. L'énergie magnétique aura varié d'une quantité 
ATe qui dépend du déplacement donné aux balais. L'observateur sera ainsi conduit à considérer 
existence de dérivées: 

d TE , dTE 


+ 0; — + 0 
do de 


(20) 


et à definir les paramétres comme des fonctions des coordonnées o. 

On doit remarquer que les points de vue de l'observateur des valeurs moyennes et de celui des 
valeurs instantanées sont opposées en ce qui concerne les relations (19) et (20). 

En fait, pour ce dernier observateur, les dérivées (20) sont nulles et les (19) différentes de zero. 

On a considéré jusqu'ici que l'armature était lisse, ce qui nous a conduit aux relations (19). Par 
contre, si larmature possede une saillance magnétique, les dérivées (19) ne seront plus nulles pour 

. a 4 - 4 7 * * 
l'observateur des valeurs movyennes. Elles joueront alors le róle des dérivées introduites dans 
ou 


les équations (12): 
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do Te dTe dTe STE 


2 am aa 


Dans les proprietés qu'on vient de définir, on reconnait tout de suite celles du modele égquiva- 
lent à une machine à commutation qu'on a introduit dans une publication antérieure [3]. 

Si on reprend les équations (14) et (18) et qu'on les écrive du point de vue de Iobservateur 
des valeurs moyennes, elles prendront la forme: 


d*T'e di PTE a 9*T Nica OF 
(22) nshati o tia ti da IES 
b Oia dip dt u dia dO bu dia OPbu u ja Pau 

Y o T' 97 

(23) (1 u) em = - É + ps Cas 


d Ou a O Pau 


L'équation électrique (22) peut être écrite sous la forme plus simple: 


d odTe mau O TE 


24 Va = Rai 
i no dt dia u lia O Pau 


On retrouve ainsi les résultats qu'on avait établis directement pour le modele équivalent des ma- 
chines à commutation [3]. 

IH est à remarquer dans ces résultats que les forces électrodynamiques et les f.e.ms. 
Z = r rm ne sont pas calculées comme des forces d'inertie généralisées. 
u la au 

Ce fait montre tout de suite que les équations de Lagrange ne sont plus valables pour Iobser- 
vateur des valeurs moyennes alors qu'elles I'étaient pour lobservateur des valeurs instantanées. 

Pour trouver comment Iobservateur des valeurs moyennes pourra interpréter le comportement 
de la machine d'un point de vue dynamique, il faut se reporter aux relations (11). 

Ces relations traduisent le principe de conservation de I'énergie pour 'observateur des valeurs 
moyennes. 


En fait, considérons la relation 
(25) 2 Yaia = Te + Te 


ou ne figure pas I'énergie associée aux sections en commutation, laquelle est infinitésimale et ignó- 
rée par lobservateur des valeurs moyennes. 
Cet observateur peut établir à partir de (25), par dérivation et aprês sommations, l'équation : 


o d Va 7, Va 
26 p + Z a [O =— 
a adu d Pbu a ad o) o 
d Te dTe Te - dTe . 
=— =D sz O + Z —— 2—— 6 
b, u O Pu a + u dou E bu O pu “bu ha u dO, n 


Si on compare ce résultat avec les equations (22) et (23), pour ip==c!., on vérifiie que, en 
valeur moyenne sur intervalle de commutation, une partie des travaux électriques et mécaniques 
est compensée par la variation d'énergie magnétique du systême. Ces travaux pourront être associés 
au travail des forces d'inertie généralisées du systême. 
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Nau STE á á 
—— et encore le travail, avec la vi- 
la au 


Pourtant, il reste à compenser le travail des f.e.ms. 


tesse nau , des forces électrodyinamiques. 

La relation (11) traduit précisément le fait que ces travaux électriques et mécaniques se com- 
pensent deux à deux. 

Pour I'observateur des valeurs moyennes, ces travaux seront dus aux liaisons du systeme. 

En fait, reprennons le point de vue de lobservateur des valeurs instantanées et établissons un 
bilan énergétique de la machine sur Iintervalle de commutation. 

Admettons, comme en I'a déjà fai précédemment, qu'un circuit b est déplacé d'un angle à Opy 
égal à l'intervalle mécanique de commutation, tous les autres circuits restant fixes. Les courants sont 
maintenus constants, à I'exception du courant de commutation ipy lequel sera bien obligé de changer 
de (—ib) par suite des liaisons intérieures, même si cela arrive brutalement à la fin du déplace- 
ment 4 Op. 

Etudions la variation pendant le déplacement de la fontion 


(27) 2 iabva + 2 iau Vau = Te + T'E 
a 


a,u 


Cette fonction joue, pour Iobservateur des valeurs instantanées, le rôle de la fonction (25). 

Il faut remarquer que les termes ia Yau ont une valeur qui est négligeable pour les deux 
observateurs. Pourtant leurs dérivées ne sont pas négligeables pour l'observateur des valeurs instan- 
tanées lequel est ainsi obligé de les considérer dans l'analyse de fonctionnement de la machine. 

La variation de (27) pendant le déplacement considéré s'écrit : 


ib 
"Abby O cas ss ; 
Rr” y, A, ; 
(28) 2 ia O Ya d Opy + 2 ia Ata dipy + 
> ” d Oby ; ib a 0 iby 
2 
ip E Pet, ip 
2 — maio PS Oby , , 
2 o T'E 2 dT'Ee À 
+ Yby d iby = À Te + - dOpy + mê São di 
) à Ta d Opy ip O iby 
2 2 


lorsque Ion néglige les variations Alau infinitésimales par rapport aux autres. 
Si on considêre A Op, suffisament petit pour les dérivées figurant en (28) soient constantes, 


on peut écrire: 


E É jT'e 
(29) Eiao a A Opy + Z ia Ya (— ip) + Ypo (— ib) = A Te + Po, 4 Opy Ed (— ib) 
a d py a d iby d by 


Par simple observation du systême, on sait que À TE = O et que les variations totales du flux 
à travers les circuits électriques, à I'exception du circuit b, sont aussi nulles. 

Si Ion tient encore compte du fait que d T'g / dipy représente le flux Yby à travers les deux 
sections en commutation, l'éguation (29) pourra s'écrire : 


(30) ib A, Yb — ib Yby aus Phy a! Opy e Yby ib 
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les deux membres étant d'ailleurs nuls, parce que la variation totale de la fonction (27) est nulle 
pour le déplacement envisagé. 

L'équation (30) montre que le travail de la force électrodynamique se retrouve tout entier 
donné dans le circuit électrique b par le terme ip Ap. 

Dans ce terme, À Jp = py représente la variation totale du flux Yp, laquelle est due au dépla- 
cement A Opy et à la variation (— ib) du courant iv. 

L'observateur des valeurs moyennes qui trouve directement ces travaux électriques et mécani- 
ques, les attribuera aux liaisons du systême, qui présente ainsi pour lui la possibilité de réaliser 
une conversion directe d'énergie. 

La condition (11) représente alors pour cet observateur, une condition du travail nul des liaisons du 
systeme, qu'il écrira dans le cas général sous la forme: 


JT'E dT'e 1 
31) dO, — dg oa EE dg 
dg : o dio Ato 


:Dbv 


ayant désigné par qb la charge mésurée à travers une section transversale du circuit extérieur b, 
Cette condition est intrinsêéque au mouvement réel du systême et rêsulte du seul fait que la 
configuration de celui-ci se répête. Les coefficients de (31) prennent automatiquement la valeur ne- 
céssaire pour que la relation soit vérifiée quelles que soit les variations dOu, dopy et dgp. 
Evidemment, la situation est différente quand on considêre des déplacements virtuels. Les coef. 
ficients de (31) ayant alors une valeur déterminée, seulement certaines combinaisons 20, dºpy et dgp 
conduisent à un travail nul des liaisons intérieures. 


Si Ion introduit la valeur des coefficients calculés sur la trajectoire réelle, les conditions (31) 
deviennent 


(32) > 10, — à, dt, — (O, — pro! “q, = O 


qui est la condition utilisée dans une publication antérieure [3]. Cette condition permet de modifier 
les équations de Lagrange de façon qu'elles puissent décrire le comportement moyen des machines 
à commutation Les résultats ainsi obtenus conduisent aux équations (24) et (23). 


Iv. CONCLUSION 


L'étude qu'on vient de réaliser montre pourquoi on peut réduire une machine à commutation 
discontinue à un modêle équivalent à commutation continue. 

Ce passage au modéle équivalent simplifie beacoup l'étude du fonctionnement de la machine. 
La raison en est que, dans I'étude de la machine réelle, nous sommes obligés de tenir compte de la 
façon dont sont réalisées les liaisons internes de la machine, alors que dans l'étude du modêle équi- 
valent, nous tenions seulement compte du mode d'action de ces liaisons. 

Ce fait ne doit pas être utilisé pour refuser intérêt de l'analyse instantanée de la machine 
réelle. En fait, il y a des machines oú les aproximations introduites en considérant la période Ate 
infinitésimale et I'énergie associée aux sections en commutation négligeable, paraissent peu précises : 


cela peut arriver dans les machines modernes à commutation électronique qui possêdent en général 
un nombre três réduit de sections. 
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Dans un cas comme celui-ci, la méthode du modele équivalent donnera des renseignements 
sur le comportement moyen de la machine, ce qui est important pour prévoir les diagrammes énergé- 
tiques du systême. Mais une analyse approfondie doit considérer le comportement instantané de 
la machine. 
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TECNICA XIV 


C. D.U. 519:622,1 


PRESENTATION DES VARIABLES RÉGIONALISÉES (*) 


RESUMO 


O autor apresenta, numa 1.º parte, uma exposição sobre 
as representações transitivas das variáveis regionalizadas, 
com referência especial ao covariograma transitivo e sua 
significação estrutural. Indica-se o cálculo da variância 
de estimação associada a uma malha de amostragem, e 
fus-se uma aplicação à estimação de uma área minerali- 
zada, Numa 2.º parte, indica-se a representação das variá- 
veis regionalizadas por esquemas intrínsecos, com referên- 
eia especial ao variograma intrênseco e sua significação 
estrutural, Indica-se o cálculo das variâncias de extensão, 
e dão-se as fórmulas gerais e métodos de cálculo da variân- 
cia de estimação associada a uma rede de trabalhos de 
reconhecimento, 


I — INTRODUCTION 


Variables régionalisées 


Un grand nombre de phénomeênes naturels se 
présentent à Ihomme sous forme régionalisée : 
ils se deploient ou se distribuent dans espace. 
De tels phénomênes peuvent se caractériser, lo- 
calement, par certaines grandeurs qui varient 
dans l'espace, et constituent, par conséquent, des 
fonctions numériques (ordinaires). Ce sont de 
telles fonctions numériques que nous appelons 
des variables régionalisées : il s'agit lá d'un terme 
neutre, purement descriptif, antérieur, en parti- 
culier à toute interprétation probabiliste. 

Comme exemples simples de variables régio- 
nalisées, citon entre autres: 


— la teneur, dans un gisement minier, 

— le rendement à lhectare, dans une région 
céréaliêre, 

— la densité de population, à lintérieur d'un 
pays. 

Dans ce qui suit, j'emprunterai surtout mes 


par G. MATHERON 


Ingénieur au corps des Mi- 
nes, Docteur és-Sciences — 
Ecole Nationale Supérieure 
des Mines de Paris 


SYNOPSIS 


The author presents, in a first part, a study on the tran- 
sitive representation of the regionalized variables, with a 
special reference to the transiítive covariogramm and its 
structural meaning. A computation of the varianee of esti- 
mation associated to a camplino net is shown, with an 
exampte of the estimation of a mineralized area 

In a second part, the representation of the regionalized 
variables by probabilistic schemes 1s indicated, with a spe- 
cial reference to the antrinsec variogramm and its structu- 
ral meaning. àls0 shown 18 the calculation of the variances 
of extension, as well as the general formulae and caleula- 
tion methods of the variance of estimation associated to un 
exploration work's net, 


exemples aux sciences de la terre et à Iestima- 
tion des gisements miniers, domaines ou j'ai 
surtout travaillé, mais le champ d'application est 
évidemment plus général. 


Caractere mixte aleatoire/structure 


Le plus souvent, les variables régionalisées 
présentent un haut degré d'irrégularité. Elles ne 
sont généralement pas dérivables, ni même con- 
tinues. On ne peut les représenter graphique- 
ment que d'une maniére grossieére et approxi- 
mative, sous forme de courbes en «dents de 
scie». L'étude directe de courbes de ce genre 
serait prohibitive, et ne présenterait sans doute 
qu'un intérêt limité, en raison de leur complica- 
tion même. Il y a dans le comportement des 
variables régionalisées un aspect aléatoire, qui sug- 
gere presque irrésistiblement le recours à une 
interprétation probabiliste. 

Cet aspect aléatoire permet de comprende 
Vinsuffisance des méthodes traditionnelles d“esti- 


(*) O presente artigo é constituído pelo texto das conferências promovidas pela Secção de Tratamento de Miné- 
rios do Centro de Mineralogia e Geologia da C.E.E.N., e realizadas nos dias 15, 16 e 17 de Março de 1967 no Instituto 
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mation des gisements miniers. Les mineurs 
admettent que chaque échantillon prélevé repré- 
sente la teneur réelle de sa zone d'influence. Ils 
obtiennent bien ainsi une estimation globale uti- 
lisable, mais ne peuvent évidemment pas l'assor- 
tir d'un calcul d'erreur (d'une variance d'estima- 
tion). Il est clair qu'entre I'échantillon et sa zone 
d'influence il n'y a pas identité, mais seulement 
corrélation, d'autant meilleure que la zone d'in- 
fluence est plus petite et que la minéralisation 
est elle-même plus continue. 

De même encore les méthodes d'interpolation 
fonctionnelle surestiment, en général, de façon 
inadmissible le degré de continuité du phêno- 
meêne représenté. Par quatorze points expéri- 
mentaux, on peut toujours faire passer un poly- 
nome du treiziême degré. Mais, en général, ce 
polynome ne refléte pas le moins du monde 
Vévolution réelle du phénomeêne entre les points 
expérimentaux. 

Et cependant, sous la complication et l'irrégu- 
larité extrême que présente une régionalisation 
dans sa variation spatiale se dissimule, en géné- 
ral la structure d'un phénomêine naturel. C'est de cette 
structure spatiale qu'a leur tour ne peuvent pas 
rendre compte les méthodes purement statisti- 
ques. Quand on classe les échantillons sous 
forme d'histogramme, on fait, par lã même, 
abstraction de l'endroit ou ils ont été prélevés: 
on détruit les structures spatiales, qui constituent 
justemente laspect le plus important du phéno- 
meéne. Dans un gisement minier, par exemple, il 
ne suffit pas de savoir avec quelle fréquence se 
répete une teneur donnée. On a besoin, aussi, 
d'apprécier le degré de continuité de la minéra- 
lisation, de connaitre la taille et la position des 
zones exploitables, etc... 

Ce n'est pas lã seulement une vue de Iesprit. 
[ est três clair que lerreur que Ion commet 
lorsque Ion estime les reserves d'un gisement à 
partir d'un rese u donné de travaux de reconnais- 
sance depend au plus haut point du degré de 
continuité de la minéralisation. Ainsi, la continuité 
d'une minéralisation constitue l'un de ces carac- 
teres structuraux dont lhistogramme statistique 
est incapable de rendre compte, et qui doivent 
cependant imperativement être pris en charge par 
une theorie réaliste des variables regionalisées. 
H y a bien d'autres caractéres analogues, nous 
les mentionnerons au fur et à mesure que nous 
les rencontrerons. 
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Le langage des fonctions aleatoires 


Il fallait donc adopter un mode de formulation 
synthétique, capable de rendre compte à la fois 
de ces deux caractêres contradictoires des varia- 
bles regionalisées, de leur aspect à la fois aléa- 
toire et structuré, et permettant bien entendu 
aussi de poser correctement et de résoudre des 
problémes essentiellement pratiques, como celui 
de I'erreur que Ion commet lorsque Ion estime 
une variable regionalisée à partir d'un échan- 
tillonnage fragmentaire. 

Ce langage adequat-à la fois probabiliste et 
structuraliste — ce sera celui des fonctions aléatoires. 
Mais il se posera un problême methodologique 
assez grave, sur lequel nous reviendrons ci-des- 
sous. Ce langage des fonctions aléatoires impli- 
que, en effet, certaines hypotheses fondamenta- 
les, de nature probabiliste, et nous devrons pro- 
ceder à un examen critique de leur sens et de 
leur valeur. Autrement dit, nous devrons exa- 
miner si, en parlant de fonctions aléatoires, 
nous disons réellement quelque chose de sensé, 
ou bien si, par hasard, nous ne parlerions pas 
pour ne rien dire. 

Avant, donc, d'introduire et d'examiner ces 
hypotheses probabilistes, il est sage et utile de 
regarder au préalable jusqu'ou il est possible 
d'aller sans introduire de telles hypothêses, c'est 
à dire en utilisant uniquement des méthodes de 
nature purement geometrique (non probabiliste). 
En fait, nous allons le voir, il est possible d'al- 
ler assez loin, et, au fond, on pourrait fonder la 
geostastique en se passant completement d'hypo- 
theses probabilistes. 

Nous exposerons donc, en premier lieu, ces 
méthodes geometriques, que Ion appelle métho- 
des transitives, et nous consacrerons une deuxiême 
partie à leur correspondant probabiliste. Il se 
trouve que — par une convergence remarquable — 
les resultats ultimes auxquels conduisent ces 
deux groupes de méthodes sont pratiquement 
equivalents. C'est lã une circonstance três ras- 
surante sur le plan methodologique. Elle signi- 
fie, en effet, que la part d'arbitraire qu'entraine 
inevitablement introduction d'hypothêses proba- 
bilistes est en fait três réduite. 


IH — LES METHODES TRANSITIVES 


Nous examinerons d'abord le cas particulier, 


PECNICA N.º 365 


de nature purement geometrique, oú il s'agit 
simplement de décrire de maniêre adéquate la 
forme d'une surface ou d'un volume. Nous intro- 
duirons ensuite la notion générale de covario- 
gramme transitif, et nous examinerons sa significa- 
tion structurale. Nous montrerons ensuite que 
cet outil minimal permet à lui seul de résoudre 
le problême de l'estimation. 


Le covariogramme geometrique 


Considérons le phenomêne de transition le 
plus simple que Ion puisse imaginer; présence 
ou absence d'un caractêre. Soit, par exemple, 
une formation geologique S d'extension limitée. 
Un sondage foré au point x la rencontre, ou ne 
la rencontre pas. Posons: 


EGO == Isixes 
O sixesS 
Il s“agit d'un phenomene unique, pour lequel 
aucune formulation probabiliste n'est possible. 
Parler de la probabilité pour qu'un point x donné 
appartienne à S n'aurait pas grand sens. On peut 
noter que tout l'interêt se concentre sur la frontiêre 
de S. En effet, k(x) est constante à lintérieur 
comme à l'extérieur de 5, et c'est au franchisse- 
ment de cette frontiêre seulement que k (x) varie, 
passant de 1 à 0, ou de O à 1. De lã le nom de 
phenomêne de transition et de méthodes transi- 
tives. 
L'aire S de notre formation est evidemment 
donnée par lintégrale 


Sm fio dx 


La valeur de S constitue un renseignement fort 
interessant du point de vue pratique. Le plus 
souvent, c'est elle que Ion cherche à estimer, à 
partir d'un réseau de sondages. Toutefois, elle 
n'apporte encore aucune information de nature 
réellement structurale. 

On peut, en effet, définir la structure d'un en- 
semble comme le systême des relations existant 
entre les éléments ou les parties de cet ensemble. 
Ainsi, nous n'aurons d'information de nature 
structurale sur notre surface S qu'à la condition 
de faire intervenir simultanement au moins 
deux points. 

Soient donc x et x + h deux points, c'est à 
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dire le plus petit ensemble structurant que nous 
puissions imaginer. Considérons alors l'expres- 
sion k x) k (x + h). Elle vaut 1Isixetx+h 
appartiennent tous les deux à S, et O autrement. 
Mais dire que x-+h appartient à S, c'est dire 
que x appartient au translaté Sp de S dans la 
translation—h. Ainsi, on a: 


Isix eSNSn 
O six|eSNS 


k (x) k (x + h) = 


Si nous intégrons cette expression en x, nous 
obtenons une fonction de h: 


(1) Kth) =| k 6) k (x + h) dx =Mes| SS. | 


qui représente la mesure (l'aire) de I'intersection 
de S et de son translaté par —h. Cette fonction 
est symétrique, puisque les deux intersections 
S5nS.n e 5NSn se deduisent une de Iautre 
par translation. On a donc: 


e (h) =XK (—h) 
K (h) <. K (0) 


Cette fonction K (h) est le covariogramme géome- 
trique associé à S. Il donne une certaine image de 
la forme de l'ensemble S. En effet: 


1— Tout d'abord, la mesure de S est K (0), 
comme on le voit en faisant h=-=O0 dans (1). De 
même, si l'on intégre (1) en h, on obtient le carré 
S de la mesure de S: 


K (0) =sS 
JK(h) dh=S? 


2—Sil'on designe par r=-|h | et parzle 
module et la direction du vecteurh, K (h) est une 
fonction de la forme K (r; 2). Dans chaque direc- 
tion «, la fonction K (r;z) est identiquement 
nulle lorsque r devient superieur à une valeur 
limite a (x) que nous appellerons portée du cova- 
riogramme K (h). La portée représente la plus 
grande dimension de S mesurée dans la direction 
«. Pour | h | > a (2), Vintersection SNS. n est 
vide. Naturellement, la portée a (x) dépend de la 
direction «, et la maniére dont elle se modifie 
lorsqu' « varie donne une image des proprietés 
d'anisotropie de la surface 5, 


3 -—Le comportement du covariogramme au 
voisinage de lorigine présente lui aussi un grand 
interêt. Il est toujours linéaire. En effet,si à r est 
três petit, aire balayée par le petit vecteur de 
module à r et de direction x dont I'origine décrit 
le contour de S vaut: 


2 [K(0)—K(3 r;2)] 


Mais, si nous désignons par 2D» la somme des 
projections sur un axe perpendiculaire à « des 
eléments d'arc du contour de S, cette aire est 
aussi égale à 2Dsir. On a donc: 


K (0) —K (3r; 2) =D,ir 
c'est a dire: 
DE K' (o , x) e DD, 


Dans le cas oú S est convexe, Ds» est la longueur 
du contour apparent de S en projection sur l'axe 
perpendiculaire à «2, et cette importante caracte- 
ristique determine la pente à l'origine du covario- 


gramme. 


Fig. 1 


Le covariogramme transitif 


Plus généralement, au lieu d'ue fonction k(x) 
susceptible de prendre les seules valeurs O et 1, 
on peut considérer une fonction f(x) represen- 
tant une variable regiónalisée quelconque. Nous 
supposerons seulement que f(x) n'est differente 
de O que dans un domaine borné. A cette fonc- 
tion f(x), nous associerons son covariogramme g(h) 
défini par: 


(2) g (h) = [ f(x) Flx + h)dx 


Cette fonction g(h) est symétrique et vérifie 
les relations: 


| bia) 
g (h) < g(0) 


Exactement comme son homologue K (h),g'h) 
reflête certaines proprietés structurales de la ré- 
gionalisation : 
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1 — Désignons par Q= f f(x) dx Vintégrale 


de f(x) —la «quantité de métal». En intégrant 
(2) en h, on obtient: 


feman=o" 


Par contre, la valeur g(0) de g(h) en h=0 
n'est pas égale à Q, mais à lintégrale du carré 


de f(x): 
g(0) = / [GJP dx 


2— A chaque direction correspond une portée 
a (2), distance au delã de laquelle g (h) est iden- 
tiguement nul. La maniére dont la portée se mo- 
difie en fonction de « reflete les anisotropies de 
la régionalisation. 

3— Le comportement de g (h) ao voisinage de 
V'origine, de son coté, donne une image du degré 
de continuité ou de régularité de la variable ré- 
gionalisée. En effet, on a toujours: 


1 fe 19 
sl)-gh)= | [E x+h)—fo)]? dx 


Donc, tout d'abord, si f (x) est partout dérivable, 
g (h) va presenter un comportement parabolique 
au voisinage de Iorigine. 

Par contre, si f(x) est continu “par morceaux” 
— c'est à dire continue et dérivable partout sauf 
en certains points constituants un ensemble de 
mesure nulle oú elle presente des discontinuités 
de premiere espéce, ou sauts — g (h) présentera 
un comportement linéaire au voisinage de Iori- 
gine: tel était le cas du covariogramme K(h) 
associé à la variable k (x) passant de + 1ã0 
à la traversée du contour de 5. 

Plus généralment, on caracterisera un g (h) par 
son developpement limité au voisinagedeh=0: 


g (h) = g (0) + ah+tah+... + oh? 


Dans ce developpement, on distinguera deux 
parties: une partie régulitre, comprenant unique- 
ment des termes entiers pairs h'”, et une partie 
irrégulitre comprenant des termes en h” () diffé- 
rent d'un entier pair) ou éventuellement des ter- 
mes du type h” log h. 

Si la partie réguliere existe seule, g (h) est in- 
définiment dérivable en h=-0, et on peut alors 
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montrer qu'il en est de même de f (x). C'est donc 
la partie irrégulitre qui caractêrise le degré de plus 
ou moins grande continuité de f (x), et avant tout 
son terme de plus bas degré. Si donc à est I'ordre 
du terme irrégulier de plus bas degré, À est un 
parametre structural donnant une mesure de la 
régularité de la régionalisation. 


L'estimation 


Nous venons de passer rapidement en revue 
les proprietés qui font du covariogramme tran- 
sitif un outil pour Ianalyse structurale. Ce n'est 
pas le seul, loin de là, mais c'est celui qui appa- 
rait le premier. Il correspond, en effet, à la plus 
petite figure structurante que l'on puisse envisa- 
ger-celle qui est constituée de deux points seu- 
lement. 

Or — c'est une circonstance remarquable — cet 
outil structural minimal suffit, à lui tout seul, 
pour résoudre le problême de l'estimation des 
variables regionalisées: problême d'importance 
pratique considérable, 

Soit, en effet, une variable regionalisée f (x), 
et un réseau de prélevements implantés selon 
une maille reguliêre a (nous utilisons des nota- 
tions unidimensionnelles, mais la transposition 
aux espaces à 2 ou 3 dimensions est immediate). 
Si y est Iun des points de prélevements, on 
connait les valeurs numériques des f (y + pa) 
pour tous les entiers p, positifs et négatifs. 

On veut estimer l'integrale: 


Q = [ f(x) dx 
à partir de l'estimateur : 
(3) Qt) =| a o Big pa) 


Cet estimateur Q* (y) est une fonction pério- 
dique, Q (y = pa) = Q* (y), puisqu'il est manifes- 
tement indifférent de prendre comme origine y 
du réseau lun ou lautre des points de préléve- 
ment. 

Si on admet que lorigine y du réséau est 
implantée au hasard à linterieur du segment 
(pour lespace à une dimension) ou du rectangle 
(pour espace à deux dimensions) de maille a, cet 
estimateur Q* (y) devient une variable aleatoire 
de moyenne Q et de variance: 


* (a) = — [ [O* () — O]? ay 
a 
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Telle est la variance d'estimation associée à la 
maille a. Pour expliciter son expression, on rem- 
place Q*(y) par son expression (3), et on 
obtient facilement : 


(4) (a)=[alZg(pa)—[s(h)dh 


On conçoit, intuitivement, qu'une estimation 
doit être d'autant plus précise que le réseau de 
prélevements est plus serré, et que la régionali- 
sation est elle-même plus reguliêre dans sa varia- 
tion spatiale. L'influence de ces deux facteurs 
est visible sur la formule (4). En effet, la variance 
d'estimation apparait comme la difference entre 
la valeur approchée et la valeur exacte d'une 
même integrale Jg(h)dh, c'est à dire comme 
Verreur commise dans le calcul approché de cette 
integrale à partir des valeurs numériques des 
g (pa). Elle est d'autant plus faible que les valeurs 
numériques disponibles sont plus nombreuses 
(que le réseau de prelévement est plus serré) ; et 
aussi, comme le montre la theorie du calcul nu- 
mérique approché des intégrales, d'autant plus 
faible que la fonction g(h) est plus réguliere, 
donc que la variable regionalisée f(x) est elle- 
-même plus reguliêre, dans sa variation spatiale- 
Mais naturellement, la formule (4) donne une 
expression mathematique précise à cette vue in- 
tuitive des choses. 

Dans la pratique, cependant, la formule (4) 
conduirait à des calculs longs, puisqu'il sagit de 
mettre en evidence la difference de deux quan- 
tités voisines Iune de Iautre, et penibles à cal- 
culer. On est donc conduit à mettre au point 
des formules d'approximation permettant un cal- 
cul plus rapide. 

On peut distinguer deux cas extrêmes, selon 
que la maille a est petite ou grande devant la 
portée, 


1 — Si la maille a est plus grande que la portée, 
les g(pa) de la formule (4) sont tous nuls, à I'ex- 
ception du terme g(0) correspondant à p=0; 
il vient donc: 


(5) S(a)=lalg(o)— [g(h)dh 


Cette expression possede une signification aléa- 
toire évidente. La maille étant superieure à la 
portée, un prelevement au plus prend une valeur 
differente de O, et tout se passe comme si ce 
prelevement unique etait implanté au hasard 
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dans un rectangle égal au rectangle de maille 
et contenant la minéralisation. La variance de la 
teneur f(y) de ce prelevement unique est alors: 


1 já 1 
—— | [f(x)]º dx — — 
al) af 
et celle de lestimateur |a| f(y) est bien donné 
par (5). 


Q 


2 —Si la maille, au contraire, est petite vis- 
-a-vis de la portée, on est conduit à rechercher un 
developpement limité de c?(a) au voisinage de 
a==0. Sous reserve d'un terme fluctuant, ou 
«Zitterbewegung», à peu prês imprevisible, mais 
de valeur moyenne nulle, on peut montrer que 
ce developpement limité correspond, terme pour 
terme, à la partie irreguliêre du covariogramme 
(la contribution à ce developpement de la partie 
reguliere est identiquement nulle). 

Ainsi, dans l'espace à une dimension, on trouve 
qu'au terme r* () + 2n) correspond le terme: 

2 (a)=A;, io 
A, étant une constante (dependant de 4). Si L 
désigne la longueur du champ mineralisé total 
à 1 dimension, et n=-L'a le nombre des son- 
dages positifs (ou plutôt la valeur probable de 
ce nombre), on trouve donc pour ce terme en r*: 


1-+1 
ad 


La variance est en 1/nl+? et non pas en 1/n, 
comme laurait suggéré une application méca- 
nique de la statistique usuelle. 

Dans Iespace à deux dimensions, si aretas desi- 
gnent les deux cotés du rectangle de maille, avec 
aí <. as, et pour un covariogramme isotrope g (r) 
(ne dépendant que du rayon vecteur r et non de 
la direction du vecteur h), on obtient une expres- 
sion de la forme : 

1 4-> 
+ C; as aí 


2+ 1». 
(6) qº(ar,as) =B> a; 


Lorsque a1 est nul, cette variance se reduit à 
son premier terme, qui ne dépend que de a> et 
représente l'erreur que Ion commet en étendant 
au gisement entier la teneur moyenne des lignes 
de plus grande densité de sondages, lorsque cette 
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teneur moyenne de lignes est parfaitement con- 
nue. Le deuxiême terme, qui depend de ai, repre- 
sente lerreur que l'on commet en estimant la 
tenneur moyenne des lignes à partir des preleve- 
ments espacés de a4 sur ces lignes. La variance 
totale est la somme de ces deux termes, de sorte 
que la formule (6) exprime un principe de com- 
position des termes de lignes et de tranche: 
tout se passe comme si les erreurs que Ion 
commet, en estimant les lignes de plus grande 
densité à partir des prelevements ponctuels, et 
en estimant le gisement lui-même à partir des 
lignes supposées parfaitement connues, pouvaient 
être regardées comme deux variables indépen- 
dantes, au sens du calcul des probabilités, et 
comme si, par suite, les variances correspondantes 
pouvaient être purement et simplement addi- 
tionnées. 

Si S designe la surface mineralisée totale, et 
n==S'a,; az le nombre des sondages positifs, et 
si on pose pg ==ai/a) <. 1, la formule (6) peut 
aussi bien s écrire: 


1-+ A l | 1 A/Z 
o? (as,as) = e) a |» 11+ a +C» E | 
hn) 7 2 


La variance est cette fois en 1 n 1+%/2, 


Application a I|"estimation d'une surface mi- 
neralisée 


Nous allons maintenant appliquer la formule 
(6) dans le cas oúla fonction f (x) = k (x) ne peut 
prendre que les valeurs O ou 1, et represente 
ainsi une surface S; le problême posé est donc 
celui de Iestimation d'une surface mineralisée S 
à partir d'un reseau de sondages à maille rectan- 
gulaire ai as (ay <. as). Le covariogramme transitif 
K (h) associé à aire S est linéaire au voisinage 
de lorigine: 


K(h) =K (0) - [h | D, 


Ds: representant, comme nous l'avons vu, la 
demi-projection des élements d'arç du contour de 
S sur un axe perpendiculaire à la direction z du 
vecteur h. 


1 — Plaçons nous d'abord dans le cas isotrope, 
c'est a dire dans le cas ou Dz =D est à peu 
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prês independant de la direction «. La formule 
(6) est directement applicable, et donne: 


cá 1 D e Log 1 | a 
Ss nº? El i É a (= 


La variance est en 1/n*2, n designant le nom- 
bre des sondages positifs. Pour utiliser effecti- 
vement cette formule, à partir des données dis- 
ponibles experimentalement, on peut estimer S 


en attribbuant à chaque sondage son rectangle 
d'influence, ce qui donne: 


S = nai as 


et D à partir du contour de la réunion de ces 
zônes d'influence (cf. fig. 2). On comptera donc 
les mombres Ni; et Ng des elements parallélles à 
ai et az respectivement, qui constituent ce peri- 
metre, et on aura D == Ny ai = N; az (puisqu'il 
y a isotropie). Par suite, il vient: 


E .%F% (No 
(7) — = = Na + 0,061 E (Na < Ni) 


2 — En general, cependant, le contour ne sera 
pas suffisamment isotrope pour que Dx» puisse 
être regardée comme constante. Il présentera, 
par exemple, une direction principale d'allonge- 
ment. Si l'un des cotés de la maille est parallêle 
à cette direction principale, ce qui sera souvent 
le cas, la formule (7) précedente reste applicable : 
en effet, prenant cette direction comme axe des 
x, si nous multiplions les ordonnées par un mo- 
dule convenable, nous obtenons une nouvelle 
figure, isotrope celle-lã, pour laquele (7) est 
valable; mais cette transformation linéaire n'a 
modifié, ni n, ni Ny, ni Ns, ni la variance rela- 
tive c4/Sº, de sorte que (7) est valable aussi dans 
le cas de la figure anisotrope initiale. 


Exemple. Sur la figure 2 ci-jointe, aire minera- 
lisée S est estimée à 10 fois le rectangle de maille 
aí az. Elle comporte un trou (une lacune). Dans 


ar ( 
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le decompte de Ni et Ns doivent figurer aussi 
bien les elements du contour intérieur que ceux 
du contour extérieur. On lit donc sur la figure: 


soit Ny = 6 
soit Ns = 4 


2D;, = 12 ai 
2D:== 8as 
D'ou, par conséquent : 


2 E 

5) 1 36 » 
—* = er fa + 0,061 | no ad 
Sº 100 Ó 4 100 


soit un écart type relatif 75/S == 11/100, et une 
fourchette d'erreur relative de + 22º/0. 


HI — LES SCHEMAS INTRINSEQUES 


Dans la premiêre partie de cet expose, nous 
avons mentionné l'existence des problêmes me- 
thodologiques graves que pose I'interpreétation 
d'une variable regionalisge comme réalisation 
d'une fonction aleatoire. Ces difficultés sont dúes 
essentiellement à Iunicité des phenomeénes na- 
turels, et à la possibilité de Iinférence statis- 
tique. 

Reprenons, en effet, les termes du probléme. 
Une variable aléatoire ordinaire Y est definie par 
sa loi de probabilité G(y) =P(Y <y). Si on 
effectue un tirage au sort selon la loi G(y), on 
obtient une valeur numérique particuliêre, par 
exemple y == 98. Naturellement, à partir de cette 
valeur numérique unique y==98, il n'est pas 
possible de reconstituer la loi G(y). Il faut pour 
celã effectuer un nombre assez grand d'expe- 
riences indépendants donnant chacune des va- 
leurs y différentes. 

Une fonction aleatoire F est une variable aléa- 
toire vectorielle à une infinité de composantes: 
chacune de ces composantes représente la valeur 
(aléatoire) F(x) prise par F en chacun des points x 
de lespace. Si Ion effectue un tirage au sort 
selon la loi de F, on obtient pour chaque point x 
une valeur particuliêre f(x): f(x) constitue ce que 
on appelle une réalisation de F. 

Ainsi, entre la fonction aléatoire F et sa reali- 
sation f, il y a la même rapport qu'entre la va- 
riable aléatoire ordinaire Y et la valeur numérique 
particuliere y = 98 obtenue à issue d'un tirage 
au sort effectué selon la loi de Y.Si donc nous 
nous donnons une variable regionalisée f(x), 
nous ne pouvons pas dire que f(x) est une fonc- 
tion aléatoire. Une telle expression n'aurait pas 
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plus de sens que si Ion disait: le nombre 98 est 
une variable aléatoire. En termes corrects, l'hypo- 
thêse probabiliste que nous désirons introduire 
doit se formuler comme suit: «la variable régio- 
nalisée f(x) peut être considéree comme une 
réalisation d'une fonction aléatoire F». 

Or, il y a dans telle hypothese, et dans les 
termes mêmes ou nous l'avons éenoncêe, un as- 
pect presque platonicien: la seule réalité digne de 
ce nom devient alors, en effet, la fonction aléa- 
toire elle-même, définie idéalement par une loi de 
probabilité donnée sur un espace abstrait. Le 
phenomêne réel, régionalisé, dont nous étions 
partis, n'est plus considérê que comme une réa- 
lisation, c'est à dire comme un reflet lointain et 
déformé, une imitation maladroite de ce modele 
idéal. 

Malheureusement, nous n'avons pas acces di- 
rectement au ciel des Idées, et nous devons exa- 
miner soigneusement si une telle hypothese a un 
sens. Il faut au minimum, pour cela, que Vinfé- 
rence statistique soit possible. Or, il n'est jamais 
possible de reconstituer la loi de la variable aléa- 
toire ordinaire Y à partir d'une seule valeur nu- 
mérique. De la même maniére, on peut craindre 
qu'il ne soit pas possible de reconstruire la loi 
de Fã partir d'une réalisation unique f. Dans les 
Sciences de la Terre, contrairement à ce qui se 
passe en physique, ou l'on peut repéter une expe- 
rience aussi souvent que Ion veut, les phenome- 
nes auxquels nous avons affaire sont uniques. 
Et c'est justement parce que l'inference statisti- 
que n'est pas possible, en genéral, qu'il a fallu 
mettre au point les méthodes transitives. Mais, 
heureusement, il y a toute une categorie de 
cas ou linférence statistique redevient possible, 
même à partir d'une réalisation unique: c'est le 
cas des fonctions aléatoires stationnaires. 

Le cas stationnaire: Pour definir une fontion 
aléatoire F, il faut connaitre sa loi spatiale, c'est 
à dire toutes les lois: 


Gt Ba pi Res j= 
= P[F(xi)< 4,F (x3) <fa,...] 


pour tous les points d'appui xi,x2, ... Xk em 
nombre k quelcongque, mais fini. Ainsi, pour k 
points d'appui particuliers, on doit en général 
estimer au minimum k moments d'ordre 1 et 


k (+) 


- moments d'ordre 2, de sorte que le 
2 
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probléme reste indeterminé. La fonction aléa- 
toire F est alors dite stationnaire si sa loi spa- 
tiale est invariante par translation, c'est à dire 
si on a, pour tout vecteur h: 


G (fr, fa ...; X4, X3 +.) — 
= G(f,fe.;xa th x + bh...) 


Cette notion de stationnarité a le sens d'une 
sorte d'homogenéité statistique dans lespace: 
quel que soit le lieu ou Ion se place, le phéno- 
mene présente les mêmes caracteristiques moyen- 
nes. On peut dire, dans un langage plus imagé, 
que le phénomeêne se répete lui-même dans l'es- 
pace. Gráce à cette répétition, on conçoit que 
l'inférence statistique redevienne possible: du 
fait que la loi spatiale est invariante par transla- 
tion, le nombre des parameétres dont elle depend 
est, en effet, substantiellement réduit. 

Nous ferons donc l'hypothêse que notre varia- 
ble régionalisée peut être considérée comme une 
réalisation d'une fonction aléatoire stationnaire : 
cette hypothese n'est d'ailleurs pas toujours 
acceptable; c'est ainsi qu'un phenomêne mani- 
festant une decroissance zonée à partir d'un 
coeur riche ne peut pas être considéré comme 
stationnaire, et me peut être décrit que par les 
méthodes transitives. Il y a cependant de três 
nombreux cas ot cette hypothese stationnaire 
est parfaitement admissible. 


Variances a priori infinies 


Dans la théorie des fonctions aléatoires station- 
naires d'ordre deux, I'outil de travail est la fonction 


K (h) =E [f (x). f(x + h)] 


qui représente la covariance des valeurs prises 
par la fonction aléatoire f(x) en deux points dis- 
tants de h. Cette covariance n'existe que si la 
variance a priori existe €galement, c'est-a-dire si 
on a 

K (0) <'co 


Or, d'assez nombreux phénomênes présentent 
une capacité de dispersion illimitée et ne peuvent 
pas etre décrits correctement si on leur attribue 
une variance a priori finie. 

Ici d'ailleurs la nature nous tend une sorte de 
piege. Lorsque Ion préléve des échantillons v 
dans un champ géométrique V, on obtient un 
histogramme à partir duquel on peut toujours 
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calculer mnumériquement une variance qui prend 
ainsi une valeur parfaitement definie. Mais cette 
variance est, en réalité, une fonction qc? (v | V) 
du support v et du champ V. Elle augmente, en 
particulier, lorsque le champ V augmente, Si les 
échantillons de taille v possêdent une variance a 
priori, celle-ci doit apparaitre comme la limite 
pour Vinfini dela variance expérimentale 2º (v | V). 
IH y a une dizaine d'années, à propos du grand 
gisement d'or du Rand, des auteures d'Afrique 
du Sud (D. G. Kricr, H.S. SicHet, De Wiys, etc...), 
ont calculé, — à partir de centaine de milliers 
d'échantillons — la variance de ces échantillons 
dans des panneaux de plus en plus grands, puis 
dans une concession entiére, puis dans le gise- 
ment du Rand dans son ensemble. Ils ont ainsi 
obtenu, expérimentalement, une relation du type: 


o? (v | V) =alog — 


La croissance de la variance observêe se pour- 
suit sans défaillance selon cette loi logarithmique 
jusqu'au dernier point expérimental pour lequel 
V/v est de Iordre de plusieurs milliards. On peut 
conclure en tout certitude qu'il n'existe pas ici 
de variance a priori finie. 


Les schemas intrinseques 


Cependant, même lorsque la variance a priori 
est infinie, il arrive souvent que les accroisse- 
ments f(x + h)— f(x) conservent une variançe 
finie. On appellera schéma intrinsêque une fonc- 
tion aléatoire à accroissements stationnaires 
d'ordre deux. Etudier une fonction aléatoire en 
tant que schéma intrinséque revient à I'étudier 
par lintermédiaire de ses accroissements, c'est-à- 
dire à une constante pres (c'est ainsi d'ailleurs 
que Jon étudie le processus de Wiener - Levy ou 
les processus poissoniens). L'outil de base qui 
remplace la covariance K (h) est le demi-vario- 
gramme, défini comme l'espérance du carré des 
accroissements : 


y (h)= - E (et + h)—f 69) | 


Si la covariance existe, on a 


[oras ipa 
K (h) = 7 (2) —y(h) 
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de sorte que le variogramme et la covariance 
constituent un seul et même outil. Il nen est 
ainsi que si y (22) existe. Si au contraire y(h) 
tend vers l'infini avec | h |, la variance a priori 
K(0) est infinie et K (h) n'existe pas. Mais le 
variogramme ;y (h) reste défini et permet d'étudier 
la fonction aléatoire (ainsi, les processus classi- 
ques de Poisson ou de Wiener-Levy ont un vario- 
gramme linéaire y(h)=a | h| qui tend vers 
Vinfini). 

La varianceo? (v | V) d'un échantillon v dans un 
champ V peut s'exprimer à laide du variogram- 
me. On trouve 


(8) avi= "e | y [yr dra — 


— a bu)ç1 8) dxa 
v JvJv 


ou, si on veut: 
(9) 2 (v|V)=F(V)—F(v) 


la fonctionnelle F(V) représentant la valeur 
moyenne de ) (h) lorsque les deux extrémités du 
vecteur h décrivent le volume V, 

On voit que, si y(h) tend vers Iinfini avec h, 
F(V) tendra vers Vinfini avec V et par conséquent 
aussi la variance de v dans V. Avec un vario- 
grame du type 


7(h)=3« log |h | 


(et en supposant v et V géométriquement sem- 
blables), on obtient 


S(v | V) = log 
v 


c'est-a-dire justement la formule obtenue expéri- 
mentalement par les auteurs d'Afrique du Sud. 

De lexpression (9) résulte une relation d'addi- 
tivité. 


P(viW=(v|W+(VIV 


dont le sens intuitif est évident: la variance d'un 
échantillon v dans le gisement V est la somme 
de la variance de v dans un panneau V' et de la 
variance du panneau V' dans le gisement V. 
On peut définir également la covariance de 
deux échantillons v et ví (situés à une distance 
fixe I'un de l'autre). On obtient une formule ana- 
logue à (9), F(v) étant remplacée par la valeur 
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moyenne de /(h) lorsque les deux extrémités de h 
décrivent respectivement les volumes v et v”. Si 
y(h) tend vers Vinfini avec h, variançes et co- 
variances sont des infiniment grands équivalents 
à F(V), de sorte que le coefficient de corrélation 
expérimental tend toujours vers I'unité : ce coeffi- 
cient n'est plus un instrument de travail adéquat, 
et c'est en définitive le variogramme lui-même 
qui exprime au mieux le réseau des corrélations 
et la structure spatiale du phénomeêne : 


1) le variogramme donne un contenu précis à 
la notion traditionnelle de zone d'influence d'un 
échantillon. Sa croissance plus ou moins rapide 
reflete la maniére plus ou moins rapide dont se 
détériore l'influence d'un échantillon sur des zones 
de plus en plus lointaines du gisement. 


2) les anisotropies se manifestent par le compor- 
tement différentiel du variograme dans les diffé- 
rentes directions de l'espace. 


3) Les phénomeênes de transition se traduisent 
sur le variogramme par des paliers et des seuils 
dont lanalyse permet souvent de mettre en évi- 
dence la superposition de plusieurs structures 
d'échelles différentes. 


4) La continuité et la régularité d'une regiona- 
ion, enfin, sont três bien exprimées para le com- 
portement du variogramme au voisinage de I'ori- 
gine. Par ordre de régularité décroissante, on 
peut distinguer quatre types: 


— Comportement parabolique au voisinage de 
Vorigine: (Fig. 3) 


Fig. 3 


y (h) est deux fois dérivable em h =. 

Ce type corresponde à une variable régiona- 
lisée dérivable en moyenne quadratique, donc à 
haut degré de régularité. 


— Tangente oblique à Iorigine: (Fig. 4) 
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y lh) est continu à lorigine, mais n'est pas 
deux fois dérivable. 


Fig. 4 


Ce type correspond à une variable régionalisée 
continue mais non dérivable en moyenne qua- 
dratique, donc déjá moins réguliere. 


— Effet de pépite: (Fig. 5) 


Fig. 5 


IH se caractérise par une discontinuité à Iori- 
gine, y (h) ne tendant pas vers O avec h. Il cor- 
respond à une variable régionalisée qui n'est pas 
continue en moyenne quadratique, donc extrê- 
mement irréguliere. 


— Variogramme plat, ou effet de pépite a I'é- 
tat pur: (Fig. 6) 


Fig. 6 


Ce cas limite correspond à une variable régio- 
nalisée purement aléatoire (mesure aléatoire à 
accroissements indépendants). 


Parte irreguliere 


Lorsque y (h) est de la forme ;y (r), c'est-a-dire 
ne dépend que du rayon vecteur r = |h |, on 
peut caractériser son comportement au voisinage 
de lorigine par un développement limité de la 
forme : 

(rn) =Zaçr" + 2C, ri 
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Sur un tel développement on distingue une 
partie réguliere et une partie irréguliere. La pre- 
miére est constituée de termes de degrés entiers 
pairs. Si elle existait seule, y (r) serait indéfini- 
ment dérivable, donc aussi la variable régiona- 
lisée elle-même qui présenterait ainsi le plus 
haut degré de reégularité. C'est donc la partie 
irréguliêre, constituée de termes du type r” ou à 
est différent d'un entier pair (éventuellment aussi 
do termes logarithmiques em r*” log r) qui dé- 
finit le type d'irrégularité de la variable régiona- 
lisée, et, dans cette partie irréguliere, c'est le 
terme du plus bas degré qui joue le rôle principal. 

On est ainsi amené à définir le degré de régu- 
larité d'une variable régionalisée comme Iordre 
à du terme irrégulier r” du plus bas degré de son 
variogramme. 


Regularisation et montée 


Lorsque Ion remplace une variable régiona- 
lisée f(x) par la valeur moyenne de f(x) dans 
un échantillon v prélevé au point x (ou par 
une moyenne pondérée plus générale) on obtient 
une nouvelle variable régionalisée dont le vario- 
gramme se déduit du y(h) initial par des opéra- 
tions de convolution, sur lesquelles nous n'insis- 
terons pas ici, mais dont l'effet est toujours ré- 
gularisant. Un exemple particulierement instruc- 
tif, et important dans les applications, de cette 
régularisation est la montée: la montée, opération 
permeitant de passer d'une variable régionalisée 
définie dans l'espace à mn dimensions à des va- 
riables régionalisées définies dans les espaces à 
n-1, n-2... dimensions, est la transposition de la 
technique miniére qui consiste à forer des sonda- 
ges verticaux, tracer des niveaux horizontaux, 
etc. ... Par exemple, si f(x,y,z) est la teneur au 
point de coordonnées (x yz) à lintérieur d'une 
formation stratiforme horizontale, un sondage 
implanté au point de coordonnées x,y de la sur- 
face topographique contient la teneur (exprimée 
en quantité de métal au métre carré) 


[ E(xyz) dz 


Cette intégrale représente une nouvelle varia- 
ble régionalisée, définie dans I'espace à 2 dimen- 
sions (la surface topographique), déduite de l'an- 
cienne par une montée d'ordre 1. 

La régularisation d'une variable régionalisée à 
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la montée se manifeste por effect suivant: si le 
terme irrégulier de plus bas degré du vario- 
gramme de la variable initial est en r*, celui de 
la variable qui sen déduit par montée d'ordre 1 
est en rH*!. Le degré de régularité augmente 
d'une unité (lorsque À est un entier impair on 
a la séquence logarithmique: log r > r — 
rº log r +... ou alternent les termes impairs 1 
et les termes logarithmiques r%* log r). 


Variance d'estimation 


Montrons maintenant comment le variogramme 
permet de résoudre un problême essentiellement 
pratique, comme celui de lestimation d'une va- 
riable régionalisée à partir d'un échantillonnage 
fragmentaire. Soit par exemple à estimer la teneur 
moyenne Z d'un domaine V: 


Eb) di 
v 


à partir de prélêvements (ici supposés ponctuels) 
effectués en n points x4...xn. On forme [esti- 


mateur 
1 
Y = —— f 
E > (xi) 


et, pour apprécier lordre de grandeur de I'erreur 
possible, on introduit la variance D* (Y — Z) de 
ja différence Y — Z. Si les variances a priori exis- 
tent, on doit avoir: 


(10) D(Y—-Z)=0z + cy —20zy 

Par contre, si 7 (h) devient infini avec h, ces 
variances et cette covariance a priori n'existent 
plus. Mais la variance d'estimation Dº (Y — Z) 
reste définie, car la différence Y — Z est une 
combinaison linéaire d'accroissements de f (x). 
On montre que cette Variance est: 


(11) Dr -2=+ 5) Y (x— x) dx — 
naV Edy 


la. a ; , 
p= e fo | 10x axar 


Li] 


Cette formule a la même signification que (10), 
A la covariance 7yz correspond le terme mixte, 
ou figurent une intégrale et une sommation dis- 
crete, tandis que q?y et c*z ont pour homologues 


am 


la somme double et l'intégrale double. Mais (11) 
est plus générale que (10), et subsiste même si 
les variances a priori n existent pas. 

On notera la structure remarquable de la for- 
mule (11), oú alternent des expressions exactes 
et approchées des mêmes intégrales. Le deuxiéme 
terme apparait comme une valeur approchée du 
premier, lui-même approximation du dernier. La 
théorie du calcul numérique approché des inté- 
grales nous fait donc pressentir que la variance 
d'estimation sera d'autant plus faible: 


— que le réseau de prélévements sera plus 
serré et plus représentatif de la géométrie 
du domaine V à estimer, 

— que le variogramme y (h) sera lui plus régu- 
lier analytiquement, et donc que la variable 
régionalisée sera elle-même plus réguliêre et 
continue dans sa variation spatiale. 


Ces conclusions sont bien conformes à ce que 
suggêre lintuition. Toutefois, la formule (11) 
conduirait à des calculs assez longs, lorsque les 
prélêvements sont un peu nombreux, et on est 
amené à utiliser des principes d'approximation: 

1) Dans Iespace à une seule dimension, et pour 
des prélêvements à maille rêguliêre a, lestimation 
du champ de longueur L==na constitué de la 
juxtaposition des zones d'influence de longueur a 
au centre desquelles sont implantés chacun des n 
échantillons, admet la variance 


DIY D=6k 


n 
cp”, variance d'extension élémentaire, est fonction 
de a et du variogramme ; (h). Si rá est le terme 
irrégulier de plus bas degré, la partie principale 
de la variance d'extension est proportionnelle à a”: 
(0) a? 


= 
E” Ss 
mai 


, On voit que la variance d'es- 
n 


timation est de la forme: 


Comme a=- 


Dit— == 


n ir 


| x 
Elle est en — et non plus en — 
nÃ + À n 


dans la statistique ordinaire. Le cas |=0, qui 
correspond a un effet de pépite, redonne bien la 


comme 


1 : 
formule usuelle en ——; cela est naturel, puisque 
n 
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leffet de pépite à létat pur correspond à des 
echantillons indépendants. Mais pour une valeur 


- M T 
usuelle comme 2==1, la variance est en Rr) et 
TL 


décroit beaucoup plus vite: dans les applications, 
ce phenomêéne permet des économies substan- 
tielles. 

2) Dans espace à plusieurs dimensions, on appli- 
que un principe d'approximation connu sous le 
nom de composition des termes de tranches et 
des termes de lignes. Exposons-le sur un exemple. 
Soit un gisement filonien reconnu par traçages 
(galeries horizontales tracées dans le plan du filon) 
à des niveaux équidistants. Chaque traçage est 
lui-même échantillonné par des prélêvements 
équidistants (on suppose que la maille a de ces 
préleévements est inférieure à la relevée h entre 
traçages). L'erreur d'estimation apparait comme 
la somme de deux erreurs : 


— celle que Ion commet en étendant aux tra- 
çages la teneur moyenne des prélêvements, 

— celle que Ion commet en étendant au gise- 
ment lui-même la teneur des traçages sup- 
posée parfaitement connue. 


“A * 1 E 
La premiére erreur admet une variance — cy (a), 
2 


ou n est le nombre total de prélévements et cjº (a) 
la variance d'extension d'un échantillon dans sa 
zone d'influence, à une dimension, de longueur 
êgale à la maille a. On la calcule comme ci-dessus- 


o e 
La deuxiême erreur admet la variance = cs (h), 


N étant le nombre des niveaux tracés, et 0%» (h) 
la variance d'extension d'un traçage dans sa 
tranche d'influence (7*» (h) se calcule à partir du 
variogramme déduit de y (h) par montée d'ordre 1). 

Si a est inférieur à h, on peut montrer que ces 
deux erreurs sont à peu pres indépandantes, de 
sorte que la variance d'estimation est: 


Ts 1 
o (a) + - 
n N 


D'(Y — Z) = 52º (h) 


En genéral, les prelevements reviennent moins 
chers que le metre de galerie. Néanmoins, il est 
inutile de les multiplier indéfiniment ; la variance 
d'estimation ne descendra pas en dessous de 


59º (h), limite assez vite atteinte dês que n 
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est grand. Il arrive un moment oú on ne peut 
plus améliorer lestimation sans travaux miniers 
supplémentaires. 

Connaissant le prix de revient du prélévement 
et du mêétre de traçage, et l'expression de 742 (a) et 
os! (h) en fonction de a et de h, un calcul 
d'optimisation élémentaire permet de déterminer 
le nombre d'échantillons que Ion doit prélever 
au meétre de traçage pour minimiser la variance 
d'estimation à dépense donnée. 

D'une maniêre gênérale, la théorie des varia- 
bles régionalisées, en donnant, sous la forme d'une 
variance d'estimation, une mesure précise de 
information disponible, fournit une base solide 
pour la recherche d'optima économiques en ma- 
tiere de reconnaissance miniére. Je ne peux ma- 
lheureusement pas développer ce point ici, non 
plus que bien d'autres applications intéressantes, 
telles que: le krigeage (ou recherche de Iesti- 
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mateur optimal), ia simulation numérique de 
gisements miniers, la théorie des milieux poreux, 
etc... pour lesquelles on me pardonnera de ren- 
voyer à la bibliographie sommaire suivante: 
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ELEMENTOS SOBRE A PRODUÇÃO E O CONSUMO DE ENERGIA 
NA REDE ELÉCTRICA NACIONAL **) 


I — Breve nota mensal 


As afluências ao conjunto do sistema, exceptuados 
os aproveitamentos do Douro Internacional que se 
mantiveram sensivelmente estacionárias (na ordem 
dos 6GWh/dia) até cerca do dia 20, decresceram, na 
parte final do mês, para um nível de apenas 4GWh/dia, 
valor da mesma ordem dos verificados no começo das 
estiagens. No que se refere aos aproveitamentos do 
Douro Internacional, registase que o caudal natural 
(Douro + Esla), depois de apresentar uma pequena 
subida em meados do mês, manteve uma tendência 
decrescente situando-se no final de Abril em valores 
da ordem dos 140 m'/s, bastante reduzidos para a 
época. 

No conjunto do sistema, a afluência do mês de Abril 
totalizou 448 GWh, correspondendo- “lhe uma probabi- 
lidade de ser excedida da ordem dos 82 0/, e um coe- 
ficiente de produtibilidade de 0,63. 


II — Elementos gerais (GWh) 


| a) Mensais Variação 

| 1966 1967, 0/, 
Produção hidráulica (Ph)... 41 6 
Produção térmica (P+)..... 0, 0 O 0 
Produção total (PT). ..... 415,2 acácia + 6 


Energia recebida de empresas 


não pertencentes ao RNC (Er) 40 2 a — 48 
Exportações (Ex) . +... ...| 00 00) 0 
Importações (1) .. ....... | 0,0) 0,0) 0 
Saldo importador (51). : 0,0 0,0 O 
Consumo em bombagem (c avo 0,0 00|— O 
Produção para con- | (1) 

sumos perman. (Pcp). 344,1 967,8 + 6,9 
Produção “para con- | | 

sumos não perman. (Pcnp).. 75,1 73,2 + 3 
Total wa eda asipise 19, 5) 441,04 5 
Coeficiante de hidraalicidade 1,66, 0,65 — 


NOTA 


(1) O aumento percentual da produção para consumos permanen- 
tes, tendo em conta a incidência dos domingos e dias especiais, é 
respectivamente de s3 e 11,4 9/9: 


HI — Diagramas de carga dos dias característicos 


4.º feira: 


1966 1967 
Produção hidráulica (Ph) MWh| 14859 | 15577 
Produção térmica ( Pr) MWh 0 O 
Produção total (PT) MWh | 14859 15577 
Trocas com [ Export. (Ex) MWh 0 O 
Espanha | Import. (1) MWh 0 0 
| Consum em bomb. hidroel. (Cb) MWh O 0 
Prod, para cons. perm. (Pp) MWh| 12279 13056 
Prod. para cons. não perm, (PonP) MWh! 2580 2521 
- “TOTAL Pi (I-Ex) MWh| 14859 | 15577 
Potência máx. MW 810 | 865 
a Potência min, MW 404 84 
” I-Ex 
g E EVEN Utiliz. da ponta horas 18,8 18,0 
TE Factor de carga 0,76 0,75 
tz Potência máx. MW/ TOL | 7,6 - 
: z p Potência min. MW 299 200 
Um E Utiliz. da ponta horas 17,5 17,3 
Factor de carga 0,73 


0,72 


ABRIL 


IV — Energia armazenada nas principais albufeiras 


— bu mês 
Albufeiras ; 

cwh | mM jo (8) 

Alto Rabagão ...... | 912,2 o 99,9 
PMS soa missa 221,5 = 99,6 
Venda Nova . » . cc. «| 191,8 JB,2 
BOMBAS Pg) ksa us ú 26,6 | 964 
COBRA a mao é vó Cos 5 32,1 97,0 
VOO puras ad 80,2 86,6 
6715 RD DR 90,3 
Castelo do Bode. . .,.. 148,4 91,0 
Guilhofrei . ..... gual 43 | 92,8 
Lagoa Comprida . .... 22,9 (2) 05,1 
Canti LARA . cousa é é 44,8 12,1 
CE dg sc esmas a 10,1 73,8 
jo ERROR E Er 9,6(3) 66,8 
Tot | com À, Rabagão . . .| 2004,4 94,8 
sem A. Rabagão . . .| 103292 90,4 


NOTAS 


(1) Coeficiente de enchimento em energia das albufeiras. 


(2) Inclui 3,4 GWh armazenados em Vale do Rossim no início do 
mês e 3,4 GWh no fim do mês. 


(3) Inclui 2,7 GWh armazenados no açude do Poio no início do mês 
e 25 GWh no fim do més. 


(*) Elementos extraídos das estatisticas mensais do Repar- 
tidor Nacional de Cargas (R. N.C,). As produções é 
os consumos dus empresas do R. N. O. representam 
cerca de 949/, dos totais do Pais, 
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RECTA DE COULOMB MAIS PROVÁVEL OBTIDA 
POR ENSAIOS TRIAXIAIS DE ROTURA 


RESUMO 


Descreve-se um método de cálculo, baseado no método dos 
minimos quadrados, para a determinação da recta de Cou- 
lomb mais provável obtida a partir dos resultados de en- 
saios trixiais de rotura efectuados sobre cilindros de um 
mesmo material, Apresenta-se um exemplo numérico 


Um dos processos, de obter em laboratório a 
recta de Coulomb de um material coerente, con- 
siste em realizar ensaios triaxiais de rotura so- 
bre provetes cilíndricos. Cada ensaio triaxial, 
realizado a uma pressão lateral constante, consta 
da aplicação crescente duma tensão axial de com- 
pressão até à rotura franca do provete. 

O Laboratório Nacional de Engenharia Civil 
(LNEC) possui câmaras metálicas para a reali- 
zação de ensaios triaxiais, que permitem a apli- 
cação de pressões laterais até 1000 kg cm”? e 
de tensões axiais até 2000 kg cm”? em provetes 
cilíndricos com 5 cm de diâmetro e 12 cm de 
altura. 

A envolvente de rotura, obtida a partir dos 
círculos de Mohr correspondentes aos ensaios 
triaxiais, poderá ser assemelhada, na maior parte 
dos materiais coerentes, a um segmento de recta 
na parte das tensões normais compreendidas 
entre zero e três a quatro vezes o valor da tensão 
normal de rotura obtido em ensaios uniaxiais. 
Neste caso a envolvente de rotura recebe o nome 
de recta de Coulomb. 

Dentro desta hipótese, os círculos de Mohr, 
quando em número superior a dois, não admi- 
tirão, em geral, uma tangente comum, pelo que 
a recta de Coulomb poderá ser obtida mediante 
a aplicação de um método estatístico tendo sido, 
no caso presente, aplicado o método dos mínimos 
quadrados. 

O presente escrito descreve o método de cál- 
culo utilizado na determinação da recta de Cou- 
lomb mais provável a partir de ensaios triaxiais 
de rotura. 
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SYNOPSIS 


The author describes a procedure, based on the method 
of least squares, for determining the most probable Cou- 
lomb's straight-line from the results of triaxial failure 
tests obtained in eylinders of one and the same material, 
A numerical instance of the method is presented. 


Fig. 1 


Sejam, fig. 1, Wi e Tri as tensões lateral e axial 
de rotura correspondentes ao ensaio triaxial de 
ordem ie [ci] o seu círculo de Mohr de centro 7; 
e raio R,, dados por: 


Ti — Ma (1) 
2 

Hime Cri — Oii (2) 
2 


Considere-se a recta CS cômo a recta de Cou- 
lomb procurada, definida pelo ponto C de coor- 
denadas (0;c) e coeficiente angular tg «, em que: 


c— coesão do material coerente 
«— ângulo de atrito interno do material coe- 
rente, 


415 


sendo: 


=c+otga, (3) 


a sua equação cartesiana. 

O segmento de recta TS definido por %j: mede 
a distância entre a recta de Coulomb procurada 
e a tangente ao circulo de Mohr que lhe é para- 
lela. 

O método estatístico que se vai aplicar na de- 
terminação da recta de Coulomb é o que torna 
mínima a soma dos quadrados das distâncias di 
de todos os círculos de Mohr obtidos experi- 
mentalmente (método dos mínimos quadrados), 
considerados com igual peso de confiança. 

Da fig. 1 tira-se faciimente, tendo em consi- 
deração os triângulos [AOC] e [AS]: 


d;==R;i—c cos «—c, sena (4) 


Quadrando a expressão (4) e fazendo o seu 
somatório em relação a todos os n círculos de 
Mohr determinados, obtém-se 


A resolução do sistema de equações trigono- 
métricas (6) e (7) pode ser efectuada graficamente, 
mediante a determinação no plano (c,z) das coor- 
denadas do ponto de encontro, I, das funções Y 
e Z. Os valores das coordenadas c e 2 do pontol 
substituídos em (3) permitem escrever a equação 
da recta de Coulomb mais provável do material 
submetido a ensaios triaxiais de rotura. 

No exemplo numérico que se apresenta a se- 
guir, foram submetidas a ensaios triaxiais de 
rotura, várias amostras de rochas provenientes 
do local de fundação de uma barragem na Grécia. 
Estes ensaios faziam parte dos estudos realizados 
pelo LNEC para o referido local, estudos esses 
que constavam de ensaios no campo e em labo- 
ratório. 

Uma das rochas estudadas foi um siltito cujos 
ensaios triaxiais, no estado natural de humidade, 
cerca de 60º/o, deram os seguintes resultados : 


QUADRO 1 


Tensões em kg cm”* 


n ” n 
A=353/=s5(R-ccosa-—s sen q)? (5) 
=| = Tensão lateral Tensão axial 
si ra de rotura 
O mínimo da função 4 é obtido anulando as Tri 
derivadas parciais de À em relação aos parâme- O receio 
tros c e x, Assim, derivando e simplificando, vem: 2 | 0 320 
sa dê Lo n 62 500 
== = > k-nccosa-sena JS q= 
Nie E! i Fer 1=0 (6) 41 800 
40 1270 
od À | pe 
Z=—— =|—nsen2a)cd— (sena 3 Ri+ 
dx 2 i= 1 
n - Com base nas tensões inscritas no quadro I e 
+ cos 24 > ci) + (cosa 3 Rj— nas expressões (1) e (2), elaborou-se o quadro 
Ei o II destinado a obter os valores dos somatórios 
n ) : . = 
<A Mena E 2) aus (7) necessários ao estabelecimento das equações (6) 
QUADRO II 


Amostra | R, 
N.º (kg cm?) (kg em 
2 160 160 
62 225 275 
41 350 450 
MPE, 560 710 
2 1295 1595 
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624 575 807 825 


R; 7; | .; 
=8) (kg? cm—?) (kg? cm!) 
25 600 | 25 600 
61 875 75 625 
157 500 202 500 
“397 600 504 100 
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As equações (6) e (7) assumem os valores par- 
ticulares: 


Y=1295—4ccosa-—1l595sen x=0 (8) 


Z=2c'sen2a-— (1295 sena + 1595cos24) c + 
+ (642575 cos «— 403912 sen 2 a) = 0 (9) 


Atribuindo valores ao parâmetro « nas duas 
equações (8) e (9) obtêm-se para cada uma, me- 
diante a resolução, de uma equação do 1.º e do 
2.º graus respectivamente, os valores correspon- 
dentes do parâmetro c. 

No quadro III estão inscritos os valores atri- 
buídos ao parâmetro « e os valores calculados 
para o parâmetro c a partir das duas equações 


(8) e (9). 


c(kg cm?) 


QUADRO III 


Valores dos parâmetros « e c 


Valores calculados para c, kg cm-2, 
a partir da equação 


Valores atribuídos 
aa 


(9) 


40 
45 
50 
55 
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kg em?) 
600 


RECIA DE COULOME MAIS PROVAVEL 
T=43+Gtg47,5º 


— 


T 


; 
w 
z 
uy 
o 
z 
E 
o) 
1 
Eu 
z 
“ 
1100 1200 Ikgem2) 
TENSÃO NORMAL, O 
Fig 3 
Com base nos valores dos parâmetros « ec Na fig. 3 encontram-se traçados os quatro 


inscritos no quadro III, determinou-se gráfica- círculos de Mohr correspondentes aos ensaios 
mente, fig. 2, a solução do sistema de equações  triaxiais realizados e a recta de Coulomb mais 
(8) e (9), obtendo-se o resultado: provável obtida pelo método dos mínimos qua- 
drados, que se acaba de expor. 
c=43 kg cm-2 q == 47,5º 
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